Decomposition of flow and nonflow in relativistic heavy-ion collisions 
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We propose a method to separate A?7-dependent and A77- independent azimuthal correlations using 
two- and four-particle cumulants between pseudo-rapidity (77) bins in symmetric heavy-ion collisions. 
The At;- independent correlation may be dominated by harmonic flows, a global correlation resulting 
from the common collision geometry. The Ar;-dependent correlation can be identified as nonflow, 
particle correlations unrelated to the common geometry. Our method exploits the r] symmetry of 
the average harmonic flows and is "data-driven." We use the AMPT and hijing event generators to 
illustrate our method. We discuss the decomposed Ar;-independent and A77-dependent correlations 
regarding flow and nonflow in the models. 

PACS numbers: 25.75.-q, 25.75.Dw 



I. INTRODUCTION 

A long standing issue in relativistic heavy-ion collisions 
is the entangled flow and nonflow in particle angular cor- 
relation measurements 1]. Flow is a global property of 
particle distributions being anisotropic with respect to a 
common symmetric direction in the transverse plane of a 
collision event (e.g. the reaction plane or the participant 
plane). Flow is thought to be composed of at least two 
contributions, the hydrodynamic type global flow domi- 
nant at low transverse momentum (pr) [2| and the path- 
length dependent jet-quenching effect dominant at high 
Pt [3|. Nonflow, on the other hand, is caused by intrin- 
sic particle correlations, such as from jet-correlations and 
resonance decays H|, and is not related to the common 
symmetric direction. In general nonflow is restricted to 
few-body correlations while flow is a global correlation of 
the entire event. 

The common symmetric direction of the overlap ge- 
ometry in the configuration space is unknown a priori, 
and is estimated experimentally using particle correla- 
tions in momentum space [11,0]. The resultant anisotro_py 
measurement is therefore contaminated by nonflow 0, [3^ 
[^. There are several methods to analyze azimuthal 
anisotropy. Broadly speaking they can be categorized 
into two classes: the two-particle cumulant method and 
the multi-particle cumulant method [6|]. For example, the 
event-plane method is effectively a two-particle correla- 
tion method because a test particle is correlated with the 
event- plane reconstructed from the rest particles [1, [13] ■ 
The two-particle cumulant method is contaminated by 
nonflow, and pseudo-rapidity (77) gaps, A77, are often ap- 
plied between the particle pair to reduce nonflow con- 
tributions. However, how much nonflow still remains is 
hard to quantify. In addition, back-to-back inter-jet cor- 
relations can still contribute considerably at large A?/. 

On the other hand, nonflow contributions in the four- 
particle cumulant method are suppressed by a high or- 
der in event multiplicity and are therefore negligible in 
heavy- ion collisions pi], [l^. However, the four-particle 
cumulant method is affected by flow fluctuations differ- 
ently from the two-particle cumulant; the fluctuation ef- 



fect is negative in four-particle cumulant while positive in 
two-particle cumulant [1, 0, Q ■ The difference between 
two- and four-particle cumulants is therefore a net ef- 
fect of flow fluctuations and nonflow. So far one cannot 
separate flow and nonflow effects from the combination 
of two- and four-particle cumulants without relying on 
model assumptions. In fact, a great deal of effort has 
been invested in studying flow fluctuations and nonflow 
based on various model assumptions (T3l - [2^ . 

In this paper we propose a data-driven method to 
separate Ary-dependent and independent correlations by 
analyzing two- and four-particle cumulants of particles 
from different pseudo-rapidity bins. The Ary-dependent 
correlation can be identifled as nonflow, and the A77- 
independent correlation may be dominated by flow. We 
lay out our method in Sec. [Hi We illustrate our method 
using the AMPT event generator in Sec. IIIII We argue 
that our method is successful in separating, to a large 
extent, the flow and nonflow in AMPT based on expected 
features for flow. We employ the hijing model to render 
additional support to our method. 



II. METHOD 

When analyzing flow, one often applies 77-gaps in two- 
particle cumulant measurements. This reduces, but does 
not completely eliminate, nonflow contributions |23, f23l - 
[25} . When analyzing jet-like correlations, one often sub- 
tracts large-Ary azimuthal correlations from small-Ary 
ones to obtain the near-side intra-jet correlation [2^, [26}. 
This assumes that flow is ry- independent. In this paper 
we extend these ideas systematically, to all possible 77- 
gaps, to hopefully extract the AT^-dependent nonflow in 
entirety without any assumptions on the 77-dependcnce 
of flow. 

We propose to analyze the two- and four-particle cu- 
mulants between two 77 bins, 77^ and 77^. The azimuthal 
cumulants can be obtained from the moments 
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For four-particle cumulant, two particles are from 77^ and 
the other two from rjp. 

The two-particle cumulant is simply the event-average 
of (2) . It is defined to be the square of two-particle cu- 
mulant anisotropy @]. It is composed of three parts, 
the average flow, flow fluctuation, and nonflow. We may 
write 

V{2} = {{2)) = v(^^^)v(^^fi) + a{J]^)a{J]fi) + a'[^^) 

+S'{Vo.)S'{v^) + S{Afj) . (3) 

where 17(77) is average flow and crij]) indicates effect of 
flow fluctuation. The term cr(ria)<T{r]p) should be taken 
simply as the definition of flow fluctuation analogous to 
the fluctuation definition in a single variable quantity x 
(i.e., {x'^) = {x)"^ + cr^); It could contain product terms 
of flow and flow fluctuations. In general flow fluctuation 
can also depend on the 77-gap between the particle pair; 
we have added this part of flow fluctuation by cr'(A7y). It 
is important to note here that nonflow, generally depend- 
ing on Ary, can also have a Ary-independent component. 
For example, back-to-back inter-jet correlations at low 
Pt may not have a A77 dependence due to the sampling 
of the underlying partonic kinematics, but may gener- 
ally be ?7-dependent. When pT increases, the inter-jet 
correlation could develop A77 dependence. In Eq. ([3]) we 
have therefore included two nonflow terms, (5(A7y) for the 
Ayy-dependent nonflow and 5' {ria)5' (rjfj) for the remain- 
ing Aiy-independent component. Generally (t'(A?7) and 
(5(A?7) may also depend on rja and rj^ which we have sup- 
pressed in Eq. ([3]) for clarity. We have also suppressed 
the harmonic order n. 

The four-particle cumulant is given by 2 ((2))^ — ((4)) 
where ((4)) is the event-averaged four-particle moment. 
The four-particle cumulant anisotropy squared is defined 
to be 3 

y{4} EE V2((2))2 - ((4)) « v{r^^)v{rip) - a{in^)a{7^p) 

-a'iAr,). (4) 

It is composed of the average flow and the flow fluctua- 
tion effect. Nonflow is negligible in four-particle cumu- 
lant [HI, Note again that Eq. (U) is rather general, 
with the a terms simply indicating effects of flow fluc- 
tuations. However, flow fluctuations need to be small 
in order for the approximation in Eq. (|4]) to be valid. 
Generally speaking the fluctuation terms in Eqs. ([3]) and 
(|4]) are not necessarily the same because they come from 
cumulants of different orders; they are the same when 
fluctuations are Gaussian [6.,^]. However, for simplicity 
we have written them as if they are the same; this does 
not affect our general argument. 

We will use short-hand notations, two-particle "cumu- 
lant flow" v{2} and four-particle "cumulant flow" w{4}, 
respectively, to stand for the A?7-independent part of the 
correlation: 

V{2} = i{2}{Tjo.)vm{r,p) 



V{A} = v{4}{rj^)v{4}{ijp) 

= v{7]a)v{r]0) - a{i]a)<j{vp) . (5) 

The average flow v{ri) is only a function of 77 and is 
symmetric about mid-rapidity for symmetric heavy-ion 
collisions, viji) = v{—ri); In this case we can write wd??!). 
The same holds for (t(|?7|) and (5'(|7y|). Take the difference 
of the two- (and four-)particle cumulants for two pairs of 
77 bins, (77a, ?7^) and (77^,-77/3), where 77^ < 77^ < 0. For 
symmetric collision systems, the A77-independent terms 
in Eqs. Q and (|1]) cancel. Namely 

Ay{2} = y{2}(77„,77^)-y{2}(7;„,-77^) 

= n2}(A77i) - y{2}(A772) = Aa' + A6 , (6) 

AV{4} = V{4}{r,^,V0)-V{4}{r,o.,~r,p) 

^ n4}(A77i) - y{4}(A772) « -Aa' , (7) 

where 

Aa' = a'{Am) - 'j'{A7j2) , (8) 
A^ = 5(A7;i)-(5(A772). (9) 

and 

A7;i = Vl3 -Va, (10) 

A7;2 = -r]p - rja . (II) 

Again the A77-dependent fluctuation and nonflow terms 
a'{Ar]) and S{Ar]) may also depend on (|77q|, I77/3I) in case 
of symmetric collisions systems; these variables are sup- 
pressed in Eqs. (|8]) and ([9]). 

What remains in Ay{4} is Act' which should yield in- 
formation on it' (77). What remain in Ay{2} are Ac' and 
AS. Given information on a' from Ay{4} and assuming 
the a' in Ay{2} is similar, one should be able to deduce 
valuable information on nonflow d{Ar]). 

Note there is no assumption about the 7;-dependence of 
flow or flow fluctuation. On the other hand, our method 
relies on the symmetry of the average flow (and the other 
A77-independent terms) about mid-rapidity, and thereby 
it applies only to symmetric collision systems. 

It is important to point out that our method essen- 
tially separates A77-dependent and A77-independent cor- 
relations, not particularly flow and nonflow. Because 
flow is a single particle property, the A77-dependent cor- 
relation may be identified as nonflow. However, since 
flow is measured by particle correlations, there may be 
a A7;-dependent flow fluctuation component. The A77- 
independent correlation should be dominated by flow. 
However, nonflow correlations are possibly present in 
back-to-back pairs of particles, such as away-side intra-jet 
correlations. This part of the nonflow is likely indepen- 
dent of A77. Therefore the decomposed A77-independent 
correlation is still contaminated by away-side nonflow. 
Consequently, the A77-dependent correlation contains 
only part of the nonflow, possibly only from near-side 
correlations such as intra-jet correlations and resonance 
decays. 
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III. ILLUSTRATION WITH EVENT 
GENERATORS 

We illustrate our method usingthe AMPT (a multi- 
phase transport) event generator [23 • The AMPT model 
uses Monte Carlo Glauber initial geometry, liberates par- 
tons from the nucleons, follows a cascade evolution of 
parton-parton scatterings, and finally freezes out into 
hadrons through a coalescence mechanism. The scatter- 
ings generate collective flow among the partons and the 
final state hadrons [H, [29| . We use all default settings 
of AMPT with string melting turned on. It was shown 
that string melting was needed to reproduce the V2 mea- 
surements [2^. We use AMPT for the purpose that it 
generates events which contain collective flow, flow fluc- 
tuations due to initial geometry fluctuations, and nonflow 
such as from resonance decays and jet-correlations. It is 
not critical for our purpose whether the physics in AMPT 
is all correct or not. 

We generate total l.I x 10^ minimum-bias AMPT 
events. For illustration of our method, we focus on events 
with impact parameter 5 < 6 < II fm, corresponding to 
approximately 10-60% collision centrality. This is de- 
scribed in Sec. HiraimC] In Sec. HITdI we check our 
method against central AMPT events with < 6 < 5 fm 
(approximately top 10% centrality). We also use the Hi- 
JING event generator 3l| to check our method, hijing 
does not have flow but contain significant nonflow, and 
thus provides a critical test to our method. We gener- 
ate total 4.7 X 10^ minimum-bias hijing events. This is 
described in Sec. IIIIEl 

In both models, strong and electromagnetic decay res- 
onances are made to decay. For our study we use all par- 
ticles with pt < 2 GeV/c in our calculation of azimuthal 
cumulants. We include neutral particles for the AMPT 
study which significantly helps to increase statistics for 
the four-particle cumulant. We include only charged 
hadrons for the hijing study because it is relatively cpu- 
cheap to generate more hijing events. 

Figure [1] shows the second-harmonic four- and two- 
particle cumulant results from AMPT events of 5 < 5 < 
11 fm, V2{4} (left panel) and V2{2} (middle panel), and 
the third-harmonic two-particle cumulant, V3{2} (right 
panel), as a function of rja and rjp of the two 77 bins. 
The dip along the diagonal Ary ~ in the two-particle 
cumulants is not understood but has been also observed 
in two-particle correlations in Ref. [s^]. The A77 ^ dip 
is also observed in hijing (see Fig.ttUlin Sec. IIIIEl) . This 
suggests that its origin may reside in hijing which is used 
as the initial input to AMPT parton cascading, not from 
final state rescattering processes. 



A. A77-dependent part of flow fluctuation 

First we take the difference of the four-particle cu- 
mulants V2{4} for two pairs of 77 bins at (77a, J?^) and 
{ila, see Eq. ([7]). The contributions of average flow 



cancel and what remains is the difference Act' between 
the fluctuation effects for 77-gaps of A?7i and A772. Fig- 
ure left panel shows the difference AV2{4} plotted ver- 
sus {r]p,r]a). 

In the middle panel of Fig. [5] the difference AV2{4} 
is plotted as a function of A?72 — A?7i = — 2?7^ for var- 
ious values of rja- The difference is not very different 
from zero, suggesting that second-harmonic flow fluctu- 
ations are nearly independent of A77. Closer inspection 
seems to suggest a tendency of a decreasing AV2{4} — 
V2{4}(A77i) — V2{4}(Ar/2) with increasing A?72 — A771. 
This implies that the flow correlation between two rj bins 
is larger when the rj bins are further apart. This is coun- 
terintuitive. 

Nevertheless, we fit AV2{4} with a linear function 



Al^2{4}-c(Ary2-A?7i) 



(12) 



for each 77^ where c is a fit parameter. This functional 
form is driven by the observation in Fig. [2] and the fact 
that AV{4} must be zero when Ar/i = A?72. We plot the 
single fit parameter c in Fig.[2]right panel as a function of 
rja and find it is consistent with a constant. We therefore 
fit all the data points in Fig.[5]middle panel with a single 
function of Eq. The fit result is c = -0.07 ± 0.02 

with a xVndf = 83/77. 

The A?7-dependent part of flow fluctuation can be in- 
ferred as 



ct'(A77) =c(A77- A77„ 



(13) 



where Ar^max is chosen to be the maximum rj-gap of 9.6 
in our study. This is shown in the left panel of Fig. |4] as 
function of {r]c„rjp). Note that the parameterization of 
Eq. is obtained from the restricted region of 77^ < 
and rja < rjjs] The full region of the {ria,rjij) space is 
extended by symmetry. 

We find Ay3{4} is consistent with zero but with large 
errors. 



B. Extracting Ary-dependent nonflow correlation 

We now turn our attention to the two-particle cumu- 
lants. Take the difference of two-particle cumulants V{2} 
between two pairs of ry bins at {rja, r^p) and (?7q., —77^), see 
Eq. (jni . The contributions of average flow and the other 
Ar/-independent terms cancel and what remains should 
be the difference in nonflow between Ar/i and A?72 (and 
any difference in the A?7-dependent part of flow fluctu- 
ation). The differences AV2{2} and AV3{2} are plotted 
in Fig. [3] (left panels) as a function of {rj(3,rja). Similar 
to the procedure in Sec. IIII Al the rest of this section 
describes empirical ways to parameterize this difference. 

As seen from the figure, Ay{2} appears to decrease 
exponentially along the off-diagonal rjp — rja and vary 
modestly along the diagonal rja+rjp. We thus slice Ay{2} 
in bins of {r]a -|- ??/3)/2, and plot it as a function of 77^ — r]a 
in Fig. [3] (middle panels). Indeed the main feature is 



FIG. 1: The square root of the second-harmonic four-particle cumulant, V2{4} (left), the second-harmonic two-particle cumulant, 
l^2{2} (middle), and the third-harmonic two-particle cumulant, l^3{2} (right) in 200 GeV Au+Au collisions with 5 < fe < 11 fm 
(approximately 10-60% centrality) generated by AMPT. 
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FIG. 2: (Color online) Left: Difference AV2{4} between two pairs at rj bins of (r]a,r]fi) and (ria,—rii3). Middle: AV2{4} as a 
function of Ary2 — A?7i = —277/3 for various values of rja- The dashed lines are linear fit to c(Ar72 — A771) for each fixed value of 
rja- Right: The fit parameter c as a function of rja. The dashed curve is a straight-line fit. A global 2D-fit of c(Ar;2 — Ar;i) 
is applied to all data points in the middle panel resulting in the solid line. The fit parameter is c = —0.07 ± 0.02 (indicated 
by the solid line in the right panel) with x'^/ndf = 83/77. Data are from AMPT simulation of 200 GeV Au-|-Au collisions with 
5 < 6 < 11 f m (approximately 10-60% centrality). 



= A771 (with a downturn at 
small A771 which is related to the diagonal dip observed in 
Fig.[T|). However, this is not the case for all values of (77^-1- 
77^)/2 = — A772/2. This is because Ay{2} is the difference 
of two-particle cumulant, must vanish at A771 = A772, and 
therefore ought to have another term corresponding to 
that at A772. This term may be negligible for A772 3> A772 
(as is the case when rja and ry^ are both far removed from 
mid-rapidity; note rja < rjfj < 0), but comparable to the 
A771 term when A772 A771 (as is the case when rja and 
Tjp are both near mid-rapidity). This explains why the 
data points for [rja + rjp)/2 are quite different from 
others where (jja +rjp)/2 is far from zero. Even for those 
other data points, the large Aryi tails (i.e. large negative 
Tja and rjp ^ and thus A771 ~ A772) of AV2 seem to 
(and should) deviate from the simple exponential drop. 



In order to gain more insights, we fit each slice of data 



points in Fig. [3] middle panels by function 

AV{2} = (a-hfeA77i)e-^"i/'^-(a + 6A772)e-'^"^/'^, (14) 

where Arji = rjp — rja is treated as variable and A772 = 
^{rjp+rja) is fixed for each slice. The linear component is 
to mimic the down-turn at small Arji . We have kept the 
two, presumably nonflow, terms of identical functional 
form, but evaluated at Arji and A772, respectively. The 
individual fit result is superimposed as the dashed curve 
in Fig. [3] middle panels for each slice of data points. 

We find that the fit parameters a, b, and d are not 
constant but vary modestly with rja +77^. This is not 
surprising because, although S(Arj) may be only a func- 
tion of A77, its magnitude likely depends on the average 
(77) of the particle pair. As an example we show the fit 
parameter c? as a function of (7^1) = (77^ -t-77^)/2 in Fig. [3] 
right panels. We find the parameters a, 6, and d seem to 
all follow a cosine curve. We thus fit them to functions 
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FIG. 3: (Color online) Left panels: difference AV2{2} (upper) and AV3{2} (lower) between (77a, 77/3) and (rja, —rjp) as a function 
of {r]p,T]a). Middle panels: difference plotted as a function of r]p — rja = Ar;i for different values of {rja + r{fi)l'i = —^ri2/'2. 
Some of the data points are scaled by succussive factors of 1.2 (indicated in the legends) for clarity. The dashed curves are 
empirical fit by Eq. p4|) for each fixed (77^ +77/3)72, while the solid curves are results of a global 2D-fit to AV{2} vs (A771, A772) 
by Eq. (|16p . Right panels: parameter d from Eq. H14p fit separately for different (77^ +r]p)/2 values plotted versus (77^ +77/3)/2. 
The results motivate the Eq. (|15p fit in the dashed curve. The d parameter from the global 2D-fit is shown in the solid curve. 
Data are from AMPT simulation of 200 GeV Au+Au collisions with 5 < & < 11 fm (approximately 10-60% centrality). 



of 

a((?7i)) = ao + ai cos[a2((r?i) - as)] , 

Kim)) = ^0 + fei cos[62((f?i) - ^3)] , 

d{{T]i)) = do + rficos[d2((??i) -^3)]- (15) 

The fit result to d is superimposed in Fig. [3] right panels 
as the dashed curve. 

Given the above observations, we now fit the entire set 
of data points in each middle panel of Fig. [3] simultane- 
ously by a single two-dimensional function 

Ay{2} = [a((77i)) + 6((r7i))A77i]e-^''i/'^«''i»- 

Him)) + 6((7;2))A772]e-^''^/'^«''^» , (16) 

treating 77^ and rja as the two variables. We call this 
"global 2D-fit." Note we have kept the two terms of iden- 
tical functional form, corresponding to the A?7-dependent 
correlations of the two pairs. The parameters are given 
by the function form of Eq. but we have written 

them as a function of the average (77) of the respective 
pair. Namely, (771) is the mean 77 of the first pair (770,77^) 
with 77-gap of A771 and (772) is the mean 77 of the second 



pair (77a, -77,3) with 77-gap of A772, 

(?7l) = {Va +?7/3)/2, 

{m) = (77a - '7/3)/2 . (17) 

We perform the global 2D-fit for both AV2{2} and 
AV3{2}. The fit parameters and x^/ndf are tabulated 
in Table m The fit results are superimposed as the solid 
curves in Fig.[3]middle panels. As can be seen, the global 
2D-fit does a fairly good job to describe all the data 
points. The 2D-fit result for d is superimposed in the 
right panel and there the 2D-fit does not seem to per- 
form well which is not obvious from the middle panels of 
Fig. O However, we note that d is plotted as a function 
of rja + Va = 2(771), but the d's in the two r.h.s. terms of 
Eq. are functions of (771) and (772), respectively. The 
2D-fit is expected to be worse than the stand-alone fit to 
d in the dashed curve because of extra constraints in the 
2D-fit. 

The relatively large x^/ndf indicates that the fit func- 
tion of Eq. ([TB)) is an imperfect description of the A77- 
dependent correlations in AMPT. The imperfection of the 
fit function is amplified by the small statistical errors on 
Ay{2}. While a fit function better describing the data 
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TABLE I: Fit parameters for Ajj-dependent nonflow correlation (plus any A?7-dependent flow fluctuation effect), 5{Ari)+a'{Ari), 
and flt xVndf by Eq. (HIJ to 10-60% (5 < 6 < 11 fm) and 0-10% (0 < < 5 fm) AMPT events and to 10-60% (5 < fe < 11 fm) 
HIJING events. 





AMPT 5 < 6 < 11 fm 


AMPT < 6 < 5 fm 


HIJING 5 < b < 11 fm 


Parameter 




S3 


S2 




S2 




ao 


5.54 ± 0.02 


1.99 ±0.02 


2.1 ±0.2 


-5.1 ±0.2 


8.44 ± 0.09 


-10.5 ±0.1 


ai 


-0.82 ±0.02 


0.84 ±0.03 


-0.4 ±0.2 


6.0 ±0.2 


7.64 ± 0.09 


-10.9 ±0.1 


a2 


1.71 ±0.02 


1.32 ±0.02 


0.85 ±0.08 


0.181 ±0.003 


0.118 ±0.002 


0.0651 ± 0.0007 


as 


2.10 ±0.06 


1.72 ±0.07 


7.6 ±1.0 


31.3 ±0.6 


22.3 ± 0.3 


50.2 ±0.6 


bo 


9.7 ±0.1 


8.3 ±0.1 


5.7 ±0.2 


-44.0 ±0.6 


3.6 ±0.1 


1.0 ± 0.2 


bi 


6.9 ±0.1 


-6.7 ±0.1 


2.1 ±0.2 


71.0 ±0.8 


2.1 ±0.1 


1.2 ± 0.2 


&2 


0.652 ± 0.009 


-0.321 ±0.001 


1.2 ± 0.1 


0.0178 ± 0.0002 


0.76 ±0.01 


0.9 ±0.1 


bi 


0.46 ± 0.04 


34.5 ±0.1 


-0.8 ±0.2 


45.3 ±0.6 


8.5 ±0.2 


0.1 ± 0.4 


do 


1.094 ± 0.007 


0.84 ±0.01 


0.85 ±0.03 


5.953 ± 0.071 


2.35 ±0.02 


6.2 ±0.1 


di 


0.426 ± 0.008 


-0.22 ±0.01 


-0.07 ±0.01 


-5.30 ±0.07 


1.93 ±0.02 


-5.9 ±0.1 


d2 


0.689 ± 0.004 


0.77 ±0.03 


2.5 ±0.2 


0.12 ± 0.01 


0.272 ±0.002 


0.11 ±0.03 


ds 


3.46 ± 0.03 


-0.15 ±0.07 


-1.98 ±0.05 


-0.5 ±0.3 


11.0 ± 0.1 


0.3 ±0.7 


xVndf 


2144 / 288 


898 / 288 


306 / 288 


303 / 288 


417 / 288 


142 / 66 




FIG. 4: Parameterized Ary-dependent part of the second-harmonic flow fluctuation cr2(Ar;) (left), the second-harmonic A77- 
dependent nonflow correlation S2{Ari) (middle) and the third-harmonic Ary-dependent nonflow correlation S3{Arj) (right) in 
AMPT 200 GeV Au±Au collisions with 5 < 6 < 11 fm (approximately 10-60% centrality), plotted as a function of (?7c«,??/3). 



is desirable, Eq. ([T6|) is adequate for the purpose of il- 
lustrating our method, because the deviations of the fit 
function from the data points are small compared to the 
magnitudes of the data points themselves, as shown in 
the middle panels of Fig. [3] 

From the fitted function Eq. (ITBl) to the two-particle 
cumulant difference, we deduce the AT^-dependent corre- 
lation itself (plus any A77 dependent fiuctuation effect) 
as 

d{Ar]) + a'{Arj) = [a((ry)) ± &((77))Ary]e~^"/''(<"» . (18) 

Mathematically there is an arbitrariness in this 
procedure-an arbitrary function in {rja^rjp) is allowed 
in the deduced Ary-dependent correlations, which would 
be cancelled in the two-particle cumulant difference in 
Eq. ([TS]). Physically, however, d{Ar}) should depend on 
A77 and any Ajy-independent function only in {ria,r]p) is 
included in V{2}. Moreover, nonfiow, corresponding to 



the A?7-dependent correlations, should vanish with very 
large ?7-gap- In Eq. (jlSp we have defined zero nonflow 
at A?7 = 00, but nonflow is practically zero with rj-gap 
between the forward and backward beam rapidities. 

Figure 2] shows the parameterized second- and third- 
harmonic Ary-dependent correlations, presumably the 
nonflow ^2(^77) and Ss^Arj), as function of (riajrip), in 
the middle and right panels, respectively. Again, the pa- 
rameterization of Eq. is obtained from the restricted 
region of 770, < and rja < rjp, and is extended to the 
full (jyaj??^) space by symmetry. The variation of the 
A?7-dependent correlation amplitude along the diagonal 
is the result of the cosine variation in the flt parameter 
a{{r])). The peaks of 6 at {r]a,r]j3) = (0, 0) probably indi- 
cate excessive nonfiow correlation pairs at mid-rapidity. 
The reason for the shoulders may be due to an interplay 
between the rapidity distributions of the same-side corre- 
lation strength and the particle dN/drj of the underlying 
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background event; the (5 is a per-pair quantity, normal- 
ized by the total number of particle pairs between a given 
pair of ry bins. 

The modest variation in the fit parameter d{{ri)) shown 
in the right panels of Fig. |3] does not show up clearly in 
the presentation of Fig. The fit parameter b is to take 
into account the apparent artifact of the Arj ^ dip in 
AMPT, and is not very important to our study. 



C. Factorization of the Aiy-independent correlation 

The A?7- independent cumulants, V{2} and V{4}, 
are obtained by subtracting the parameterized Arj- 
dependent correlations from the raw two- and four- 
particle cumulants V{2} and V^{4} in Fig. [TJ As men- 
tioned in Sec. |IT1 because flow is a single particle prop- 
erty, the A?7-dependent component can be attributed to 
nonflow. Consequently the Ary-independent components 
V{2} and V{4} are dominated by flow. We note that 
V{2} and V{4} include effects of the Ary- independent flow 
fluctuation. The V{2}, in addition, may include any A77- 
independent nonflow contributions. The Ary-indcpendent 
nonflow may arise from back-to-back inter-jet correla- 
tions, however in the low px region studied in this work 
it should be small. In the following discussion, for sim- 
plicity, we refer to V{2} and V{A} as "flow" cumulants 
or "flow" and the Ary-dependent 6 as "nonflow," unless 
otherwise noted. It is important to keep in mind that 
the former (that we called flow) may contain away-side 
nonflow while the latter (that we called nonflow) may 
not contain the away-side nonflow. Our nonflow may be 
primarily caused by near-side correlations that depends 
on Arj. 

The results of V{2} and V{4} are shown in Fig. [5] In 
this section we investigate the properties of the obtained 
V{2} and V{4} to check whether they resemble what we 
would expect for flow. Figure[6]lower panels show slices of 
the obtained Ary-independent correlations in Fig. [5] at se- 
lected values for one of the r] bins (we call this rjief)- The 
data are plotted versus the value of the other 77 bin (we 
simply call rj). The data points follow a smooth pattern. 
The curves for symmetric ryrof values are consistent with 
each other as flow should be. The magnitudes of the data 
points depend on the ri-^cf value because the flow is not 
constant as a function of 77. For comparison, the slices of 
the 2D plots of the raw cumulants from Fig. [T] are shown 
in the upper panels of Fig. [S) The data are not symmet- 
ric for symmetric ?7iof values. The peaks (bumps) are due 
to nonflow contributions that are strongest at small A77. 

On the other hand, the slices of the four-particle cu- 
mulants, with and without subtracting the Ar^-dependent 
fluctuation effect, are similar as shown in the left panels 
of Fig. [SI This is because the A7y-dependent part of flow 
fluctuation is small and can be safely neglected (see Fig.|3] 
left panel). 

Next we check factorization of the flow that might be 



obtained from the diagonal elements in Fig. [S] by 



w('7rcf) = V'('7rcf,?7rcf) ^ Y I V'(?7rcf , ?7rcf ) | • (19) 

Because V{r]rcf, ?7rcf) can be negative due to fluctuations, 
we take the square root of its absolute value but retain its 
sign following a recipe from ATLAS . Note in Eq. 
we have suppressed the harmonic number subscript n and 
the cumulant ranks, such as in Vn{2} and {',i{4}. We di- 
vide the results in Fig. [5] by v{rircf) to obtain v{ri). The 
results are plotted in Fig. [7] (lower panels) as a function 
of (77rcf,?7). The results show that v{r]) is independent 
of ?7rcf , indicating that the nonflow subtracted cumulants 
are factorized as v{'i]^ef)v{r]). (The bed-headboard slices 
at large |77ref | are due to division by a small, close to zero, 
value of 'D(ryicf), with large statistical errors.) Again for 
comparison the upper panel of Fig. [7] shows the results 
obtained by the same procedure but using the raw cu- 
mulants in Fig. [T] Clearly no factorization is observed 
for the raw two-particle cumulants. On the other hand, 
the four-particle raw cumulant does show factorization; 
this is because the A77-dependent fluctuation and nonflow 
effects are negligible. 

Figure [8] shows r?ref slices of the 2D plots in Fig. [71 Note 
we have avoided the problematic bed-headboard regions 
at large Ir^refl (some effect is still noticeable for [Tyrefl = 
0.35). Figure [8] demonstrates more quantitatively the 
factorization of the A?7-independent correlations in the 
lower panels (all data points fall on a common curve) and 
the non- factorization of the raw two-particle cumulants in 
the upper panels ("flow" would depend on the reference 
particle used to measure it). 

It is important to note that flow factorization is 
not required by our method. The only require- 
ment of our method is that the average flow (and 
other A?7-independent terms) be symmetric about mid- 
rapidity, which must be true in symmetric collision 
systems. For instance, the obtained nonflow sub- 
tracted cumulants could have the functional form of 
u{ria)u{j]B) + w{rja)w(riij) which does not generally fac- 
torize ^2^. Therefore the factorization of the obtained 
A77-independent correlation is a result of the analysis us- 
ing our method, which indicates that our method may 
have succeeded to separate the major parts of flow and 
nonflow in AMPT. 

To conclude this section, we have shown that the 
ATy-independent correlations in AMPT obtained by our 
method satisfies what we would naively expect for flow- 
The diagonal element of the AT^-independent correlation 
does not depend on the reference particle rj^cf used to 
measure it. This is so despite that the AT^-independent 
correlations include flow fluctuation effects and possibly 
A77-independent nonflow contributions. The factoriza- 
tion may indicate that the flow fluctuation is proportional 
to the average flow and the nonflow effect from back-to- 
back inter-jet correlations is small compared to flow in 
the studied low-pT region. 
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FIG. 5: The second-harmonic four-particle cumulant subtracted by the parameterized Ayy-dependent flow fluctuation effect: 
V2{4} (left: panel). The two-particle cumulants subtracted by the parameterized Ar;-dependent nonflow correlation (and A77- 
dependent flow fluctuation effect): second-harmonic V2{2} (middle panel) and third-harmonic V3{2} (right panel). The results 
are plotted versus {riajTip). Data are fi-om AMPT simulation of 200 GeV Au+Au collisions with 5 < 6 < 11 fm (approximately 
10-60% centrality). 
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FIG. 6: (Color online) Slices of the 2D plots of Vn in Fig. [T] (upper panels) and those of Vn in Fig. [5] (lower panels) at selected r/ref 
values. Data are from AMPT simulation of 200 GeV Au-|-Au collisions with 5 < 6 < 11 fm (approximately 10-60% centrality). 



D. Test with AMPT central events 

We have been focusing on medium central events from 
AMPT where flow is significant. In this section we check 
the approximately 10% most central events from AMPT 
cutting on impact parameter < 5 < 5 fm. 

The statistical errors on the central data are notice- 
ably larger. The statistics in V{A} are too poor to tell 
whether there is a significant Ar^-dependent flow fluctu- 
ation effect. However, it seems reasonable to assume the 
effect is negligible based on the results from the 10-60% 



centrality. 

Given the statistical errors, the nonflow fit model of 
Eq. ((T6l) seems adequate. The global 2D-fit gives a 
X^/ndf ^ 1; see Table [H Figure |9] shows the central col- 
lision results from our method. The left panels show the 
fitted Ayy-dependent correlation, S{Ari) + (T'{Ari). 

Figure [ni middle panels show the extracted V2{'2} and 
'!33{2} as a function of rj and rj^of, similar to results in 
Fig. [T] The Ary-independent correlation results are again 
approximately independent of r/icf; the bed- headboard 
regions which show up in both W2{2} and wslS} are due 



FIG. 7: Study of factorization. The upper panels show Vn{r], rj^d) / \/Vn{i]ici, Vrci) with the raw cumulants Vn from Fig. [T] 
The lower panels show Vnir]) — Vn{ri,rjrci) /vn{rircf) with the Vn results in Fig. [S] The lower panel results are r/rcf independent, 
indicating factorization of the decomposed Ar;-independent two-particle cumulant Vn{ri,rjrci)- Data are from AMPT simulation 
of 200 GeV Au+Au collisions with 5 < 6 < 11 fm (approximately 10-60% centrality). 



to the division by the small magnitude (with large errors) 
at large jryref | in central collisions. The rirei independence 
indicates factorization of the Ayy-independent correlation 
that is expected for flow. The T^rcf independence is more 
quantitatively demonstrated in the right panels of Fig. [3] 
where the slices of 'D„{2} in ?7i-cf are depicted; all results 
fall on top of each other. 



E. Test against the hijing model 

In this section we test our method against the hi- 
jing event generator [sij. hijing is a quantum chro- 
modynamics inspired phenomenological model focusing 
on hard parton scattering processes, with particle pro- 
duction from soft processes parameterized. Since it does 
not have final state interactions, it does not generate col- 
lective flow. However, it contains significant nonfiow, 
mainly from jet-correlations. Thus hijing provides a 
critical test to our method. 

Figure [TO] shows the two-particle second- and third- 
harmonic cumulants, V2{2} and V3{2}, in the upper-left 
and lower-left panels, respectively. The main feature of 
V2{2} is the diagonal ridge along rja = Vi3j characteris- 
tic of nonfiow. In Vi3{2} the diagonal ridge seems to be 
present only in the limited mid-rapidity region. At large 



rapidities along the diagonal, V3{2} is negative. 

We follow the same procedure described in Sec. IIIIBl 
taking difference between two pairs of 77-bins, (77a, r;^) 
and {rja, —rjp). We again fit the difference Ay{2} by the 
2D function of Eq. (IT6)) . The fit parameters and x^/ndf 
are tabulated in Table HI We deduce the Ary-dependent 
nonfiow correlations from the 2D-fit. The deduced re- 
sults are depicted in the middle panels of Fig. [TOl The 
second-harmonic S2 is about an order of magnitude larger 
than the third-harmonic S3, in line with the two-particle 
cumulants in the left panels. 

Subtracting the parameterized Ajy-dependent nonflow 
correlation, we obtain the remaining Ary-independent 
correlation, 14i{2}, as a function of (770,, 77^). Because it is 
known that HIJING does not generate collective flow, the 
obtained V^„{2} are in fact nonflows due to back-to-back 
inter-jet correlations, 5'^{r]a)6'^{rifs); see Eq. We show 
V2{2} in the upper-right panel of Fig. [TO] It is peaked for 
particle pairs at mid-rapidity (0,0) and decreases towards 
large rapidities. In the lower-right panel, we show in open 
circles the projection of V2{2} to one of the T^-axes. This 
would be the average away-side nonflow of particle pairs 
with one particle at 77 and the other particle from the 
entire r/ range. 

Also superimposed in Fig. llOl lower-right panel in solid 
circles are the average parameterized nonflow values of 
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FIG. 8: (Color online) Slices of plots in Fig. [7] at selected »7rcf, showing more quantitatively factorization in the decomposed 
Ajy-independent two-particle cumulant: = V'n(?7, ??rof )/i'n(J7rof ) (lower panels) and non-factorization in the raw cumulant 
measurement V„{ri, rj,-c{)/ •y/Vii(r;rcf , ?7rof ) (upper panels). Data are from AMPT simulation of 200 GeV Au-|-Au collisions with 
5 < 6 < 11 fm (approximately 10-60% centrality). 



82- This is the A7y-dependent near-side nonflow, averaged 
over all partner particles and plotted as a function of rj 
of the test particle. The band indicates the uncertain- 
ties from the fit errors. This uncertainty also applies to 
the average V2{2} but anti-correlated. The magnitudes 
of the near-side nonflow ^2(^77) and the away-side non- 
fiow S'2{ria)S'2{ri/3) are similar. The near-side nonflow may 
largely come from intra-jet correlations, and the away- 
side may be dominated by inter-jet correlations. The 
similar magnitudes are then expected from low pT two- 
particle correlations from dijets due to momentum con- 
servation. See Sec. lIVBj for further discussion. 



The statistics of V^3{2} is poor and the two-dimensional 
lego plot does not reveal too much information. However, 
we show the projection of V3{2} in the solid triangles 
in the lower-right panel of Fig. 1101 The average near- 
side nonflow 63 is also superimposed with the shaded un- 
certainty band. Note the uncertainty band also applied 
to V3{2} but anti-correlated. The third-harmonic near- 
and away-side nonflows are significantly smaller than the 
second-harmonic nonflows. 



To conclude this subsection, the features of the sep- 
arated near- and away-side nonflows from the hijing 
model seem to all make sense. This renders additional 
support to our decomposition method. 



IV. DISCUSSION OF FLOW AND NONFLOW 
FROM THE MODELS 

In this section we discuss the A?7-dependent and in- 
dependent correlations of the models separated by our 
method. The former is the sum of the Ar^-dependent non- 
fiow S{Ari) and Ajy-dependent flow fluctuation a'{Arf). 
The latter is the sum of the average fiow, fiow fluctuation, 
and the A?7-independent part of nonflow 6' . Since a' [l\ni) 
is small as indicated by our results and 5' is presumably 
small at our low pT region, we will simply call the A?]- 
independent correlation as "flow" and the Ary-dependent 
correlation as nonflow in the following discussion unless 
otherwise noted. 



A. Comparison to AMPT model flow 

Strong interactions among constituents in relativistic 
heavy-ion collisions are expected to generate hydrody- 
namic expansion and convert the spatial anisotropy of 
the initial overlap geometry into an azimuthal anisotropy 
in the final-state particle distribution. This phenomenon 
is often called hydrodynamic flow which we refer to here 
as "real" fiow. This "real" average flow can be calculated 
from the AMPT model input because the initial geometry 
symmetry axis can be calculated from model information 
on the initial overlap geometry. In this section we com- 
pare the fiow results from our decomposition method to 
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FIG. 9: (Color online) Results for AMPT 200 GeV central Au+Au collisions with < 6 < 5 fm (approximately top 10% 
centrality). Left panels: Parameterized second-harmonic Ar;-dependent nonfiow correlation 82 (upper panel) and third- harmonic 
Ary-dependent nonfiow correlation 83 (lower panel). Middle panels: decomposed Ary-independent single-particle-level "flow" 
correlations {;2{2} (upper panel) and i33{2} (lower panel) as a function of rj (and »7rcf used to obtain it). The results are 
independent of rjicf- Right panels: Slices of plots in the middle panels at selected rjrat, showing more quantitatively factorization 
in the decomposed single-particle- level "flow" correlations V2{'2} (upper panel) and 'C3{2} (lower panel). 



the calculated "real" flow from the AMPT model to as- 
sess its performance. We note, however, that the spatial 
geometry information is inaccessible to experiment, and 
experimentally the flow axis has to be defined by final 
particle momenta. Due to particle correlations unrelated 
to flow, the experimental flow measurements are often 
contaminated by nonflow which we try to separate in 
this study. 

We project the W2{4} and V2{2} 2D plots in Fig. [7] 
lower-left and lower- middle panels over the entire |?7rcf | < 
4.8. The result is shown in Fig.[Tl]left panel by the filled 
circles and squares for second-harmonic two- and four- 
particle cumulant methods, respectively. These would be 
the two- and four-particle second-harmonic flows, largely 
devoid of nonflow, using all other particles as the ref- 
erence particles to measure them. The solid histograms 
(almost coinciding with the data points) are, respectively, 
the two- and four-particle cumulant flows obtained from 
the diagonal elements of the nonflow-subtracted cumu- 
lants by Eq. (fT9|) . The dashed histogram shows the diag- 
onal result obtained from the raw two-particle cumulant 
similar to Eq. ([T^ . The difference between the dashed 
histogram and the solid circles is the effect of the nonflow. 

As seen from Fig. [TT] the two-particle cumulant "52 {2} 



is larger than the four-particle '52 {4} and the difference is 
presumably caused by flow fluctuations (as the remaining 
away-side nonflow is expected to be small). However, 
the difference becomes diminished at \ri\ > 4. This seems 
hard to understand-it would imply that flow fluctuations 
diminish at |?7| > 4 if Eq. which assumes Gaussian 
and relatively small flow fluctuation, is valid. 

Figure [TT] middle panel shows the corresponding re- 
sults in central AMPT events of < 6 < 5 fm. As 
seen in Fig. [3 the 'D2{2} obtained from large ?7icf ref- 
erence particles have large uncertainties (bed-headboard 
effect). We have thus restricted our projection to within 
jyrcfl < 3.6 in obtaining the 'C2{2} results in the middle 
panel of Fig. 1111 The four-particle cumulant in central 
collisions is not plotted-The statistical precision is too 
poor to be useful. The other results are qualitatively 
similar between central and medium-central collisions. 



We calculate the real flow in the following way. We 
obtain the minor axis (participant plane) of the coor- 
dinate space overlap region using the transverse coordi- 
nates (transverse radius r and azimuth 0^) of all partic- 
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FIG. 10: Results for hijing 200 GeV Au+Au collisions with 5 < 6 < 11 fm (approximately 10-60% centrality). Left panels: 
two-particle cumulants V2{2} (upper) and V3{2} (lower). Middle panels: Parameterized ATy-dependent nonfiow correlations 
52 (upper) and ^3 (lower). Upper-right panel: 82 subtracted two-particle cumulant V2{2}. Lower-right panel: projections 
(averages) of Ki{2} (open symbols) and Sn (filled symbols) onto one rj axis. The bands are systematic uncertainties from the 
nonfiow fit errors; these uncertainties also apply to the average Ki{2} but anti-correlated. 



' V2{2}, histo:yV2{2}(Ti,Ti), dashed: yVj{2}(ii,Ti) 
' v^[4]. histo:^9^}(ti,n) 



V2{21, histo:yV2{2}(ii,il), dashed: yVj{2Kii,n) 



S 1| 

g \ o medium-central (V2{2}/V2{PP}-1)''^ 

• central (V2{2} /v^{PP}- 1)^'^ 
a medium-central [{vl{2] - ^^{i}) I [^\{2} + 02{4})f ^ 




FIG. 11: (Color online) Nonflow-subtracted second-harmonic two-particle cumulant fiow «2{2} (circles) and Ary-dependent 
fluctuation subtracted second-harmonic four-particle cumulant flow «2{4} (square) as a function of 77 in AMPT 200 GeV medium- 
central (5 < fe < 11 fm, left panel) and central (0 < 6 < 5 fm, middle panel) Au-f An collisions. Right panel: estimates of 
relative V2 fluctuation, 02 /v2, in AMPT. 



ipant gluons and quarks by [13] 

^PP =: iATAN2 fV sin(n(/)r), V sin{n(t>r)] + - , 

(20) 

where n = 2, 3 for the second- and third-harmonics, re- 
spectively. The function atan2(?/, x) is the angle between 
the positive x-axis and the vector given by the coordi- 
nates {x, y) on a plane; The angle is positive for ?/ > 



and negative for y < Q. We calculate the participant 
plane resolution by using the sub-event method and find 
the resolution is essentially 100%. We calculate the real 
participant-plane flow by 

wjPP} - (cosn(0- Vpp)). (21) 

We also calculate the real reaction-plane flow by using 
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FIG. 12: (Color online) Ar;-dependent nonflow subtracted third-harmonic two-particle cumulant flow {;3{2} as a function of 77 

in AMPT 200 GeV medium-central (5 < 6 < 11 fm) and central (0 < 6 < 5 fm) Au+Au collisions. Right panel: estimates of 
relative ?;3 fluctuation, cts/ws, in AMPT. 
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FIG. 13: (Color online) Left panel: V3{2}/v2{2} ratio in AMPT 200 GeV medium-central (5 < 6 < 11 fm) and central 
(0 < 6 < 5 fm) Au-|-Au collisions. Right panel: Ratio of i'n{2} in central to that in medium-central collisions. 



the known reaction plane in AMPT (V'rp = 0) by 

«„{RP}= (cosn(0-^Rp)). (22) 

The calculated t;2{PP} and W2{R.P} in medium-central 
AMPT collisions are shown in Fig. [TT] left panel as the 
black and red curves, respectively. Those in central AMPT 
collisions are shown in the Fig. [TT] middle panel. The cen- 
tral collision W2{RP} is flatter. The W2{PP} is larger than 
W2{RP} presumably due to the effect of the fluctuating 
participant plane about the reaction plane. It is interest- 
ing to note that ti2{PP} and ti2{RP} also approach each 
other at \ri\ > 4. 

The extracted '52{2} is larger than the calculated 
ti2{PP}, the true average V2 with respect to the partici- 
pant plane from AMPT. This is presumably because the 
extracted f2{2} should contain, besides effect from fluc- 
tuating participant plane about the reaction plane, addi- 
tional flow fluctuations caused by geometry eccentricity 
fluctuations. Small contribution from A?7-independcnt 
nonflow may also play a role in V2- 

It is commonly believed that the average two- and four- 
particle, (■02{2} + W2 {4} )/2^ equals to the true average 
participant plane ti2{PP} ^6]. This does not seem to be 
true from our results in Fig. 1111 It is also commonly 
believed that the four-particle ii2{4} should equal to the 



reaction plane W2{RP} Q. This does not seem to be true 
either. 

It is interesting to examine the magnitude of flow fluc- 
tuation. However, it is not clear how one can reliably 
obtain flow fluctuation from the ^{2} and ^{4} measure- 
ments; see the previous discussion. Nevertheless, we com- 
pute y'^ ^l|2}+l ^{ij ^ extracted from our method to give 

a feeling about flow fluctuations in AMPT medium-central 
collisions. This is shown as the open squares in Fig. [TT] 
right panel. Note that the small AT^-dependent flow fluc- 
tuation effects (see Fig.[l]left panel) have been already re- 
moved from the ^{2} and v{4} measurements. One may 
also take the difference between the "measured" "52 {2} 
and the calculated W2{PP} from model input as flow fluc- 
tuation, ^u|{2}/w|{PP} — 1. This is shown as the open 
and filled circles in Fig. llll right panel for medium-central 
and central collisions, respectively. Although experimen- 
tally inaccessible, it may give additional information re- 
garding the size of flow fluctuation in AMPT. We note 
that the away-side A77- independent nonflow is included in 
the decomposed ■0^{2} so the flow fluctuations are likely 
overestimated. 

Figure [12] left panel shows the third-harmonic two- 
particle cumulants W3{2} in medium-central and central 
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AMPT collisions. For the four-particle cumulant is 
consistent with zero with large statistical errors in both 
medium-central and central collisions. As seen in Fig. [7] 
and Fig. [HI the w„{2} obtained from large ?7rcf part- 
ner particles have large uncertainties (bed-headboard ef- 
fect). We have thus restricted our projection to within 
['?rcf I < 3.6 in obtaining the w„{2} results in the left panel 
of Fig. [T^ Figure [T^ left panel also shows the calculated 
true third-harmonic flow with respect to the participant 
plane by Eq. (|2T]) . The third-harmonic flow with respect 
to the reaction plane is found to be zero at both central- 
ities. 

Figure [T^] right panel shows an estimate of 
the third-harmonic flow fluctuation in AMPT by 
•\/w|{2}/t;|{PP} — 1. Similar to the second-harmonic 
one in Fig. [Til the estimated third-harmonic flow is rel- 
atively independent of 77, and seems smaller in central 
collisions than in medium-central collisions. The relative 
third-harmonic flow fluctuation is larger than the second- 
harmonic one. 

It is interesting to examine the relative third-harmonic 
to second-harmonic flow. The left panel of Fig. [T^ shows 
the ratio of z;3{2}/?;2{2} in medium-central and central 
collisions. The V3{2} /w2{2} ratio is larger at mid-rapidity 
and decreases significantly toward forward and backward 
rapidities. The ratio is larger in central collisions than in 
medium-central collisions. 

Figure [T^ right panel shows the ratio of v{2} in cen- 
tral collisions to that in medium-central collisions. It is 
interesting to note that the ratios are approximately in- 
dependent of 77; The centrality variation of flow does not 
depend on rj. The '53{2} central to medium-central ratio 
is close to unity-the V3 magnitude has weak centrality 
dependence, consistent with its fluctuation nature. 



B. Discussion of nonflow from the models 

In this section we discuss the extracted A7y-dependent 
nonflow S{Ari) in the context of the models. This part of 
nonflow may consist primarily of effects from resonance 
decays and near-side, intra-jet correlations in the models. 
The models do not include quantum interference effects. 

Experimentally one often applies a minimum 77-gap to 
reduce nonflow. One difficult question has been how 
much nonflow still remains after a ?7-gap cut. We may 
now address this question, keeping in mind that there 
could be a A7y-independent nonflow contribution from 
back-to-back correlations that may not be cut out by 
a 7y-gap. Figure [14] shows the average nonflow with rj- 
gap between Arjcut (the horizontal axis variable) and a 
maximum ATymax = 5- We apply a maximum 77 cut be- 
cause experimentally one is limited by detector accep- 
tance. With larger ATy^ax cuts, the results are quanti- 
tatively similar. The left and middle panels show the 
average nonflow in 10-60% and 0-10% AMPT events, re- 
spectively. The second-harmonic nonflow is larger than 
the third-harmonic nonflow. The right panel shows the 



10-60% HIJING events; the third-harmonic 63 is consis- 
tent with zero (see also Sec. IIII E)) . As expected the non- 
flow decreases with rj-gap size because of a continuously 
decreasing contribution from the near-side small Acf) cor- 
relations. However, some finite nonflow is still present 
with relatively large ?7-gap. 

Nonflow is weaker in more central collisions. This is 
partially due to the large event multiplicity because non- 
flow is normalized by the total number of particle pairs 
in an event. The reduction in nonflow from medium 
central to central collisions is approximately a factor of 
1.8. The increase in the number of participants A^part 
from medium central (Apart 130) to central collisions 
(A^part ~ 325) is about a factor of 2.5. So nonflow from 
AMPT seems to scale as l/N^l^^. This may indicate that 
nonflow in AMPT comes from hard-scattering sources, the 

magnitude of which scales with A^part ^nd is then normal- 
ized by the total number of particle pairs. 

It is interesting to note that HIJING 62 is smaller than 
AMPT 62 by a factor of 5. This at first glance seems 
counter-intuitive (if 62 is indeed associated with nonflow 
as we argue). The jet-correlation effect between low-px 
particles should be either similar or even stronger in Hi- 
JING than AMPT because final state scatterings in AMPT 
would likely destroy correlations. One would therefore 
naively expect, instead, a larger nonflow in hijing than 
in AMPT. One possible reason for the observation is 
that the parton interactions in AMPT have moved jet- 
correlated particles to lower px, so that the correlation 
strength at low pT is relatively stronger. Experimentally 
it is found that the away-sidejet-medium interactions 
soften the correlated particles [Sj] while the near-side jet 
was biased toward surface emission by the high-pr trig- 
ger particles so little modification was observed. Since 
the particles used in our study are at low px, effects of 
jet-medium interactions should be also visible on the near 
side, which is probably what is revealed in the nonflow 
difference between AMPT and hijing in Fig. [Hj 

It is also interesting to examine the magnitude of non- 
flow relative to flow. Figure [121 left and middle pan- 
els show the ratios of average ((5„) to (■u^{2}) in AMPT 
medium-central and central collisions as a function of 
the minimum ry-gap, Arycut- A maximum Arymax = 5 
is applied as before and the average (w^{2}) is calcu- 
lated in the same way as the average {6n)- The rel- 
ative third-harmonic nonflow to flow ratio {S) / {v'^) is 
a factor of 4 larger than the second-harmonic ratio in 
medium-central collisions. This is mostly due to the 
much larger V2 than in these collisions. In central colli- 
sions, {S)/{v'^) becomes comparable between the second- 
and third-harmonics. 

The right panel of Fig. [TSl shows the ratio of nonflow to 
"flow" from the hijing model. As discussed in Sec. IIII El 
hijing does not have flow; the extracted flow reflects the 
Ayy-independent away-side nonflow contributions. Thus 
the ratio reflects the relative near-side to away-side non- 
flow strengths. With a minimum ?y-gap of zero, the ra- 
tio of the second-harmonic nonflows is about unity- this 
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FIG. 14: (Color online) Average magnitude of ATy-dependent nonflow as a function of the minimum 7)-gap size (with a maximum 
r;-gap of 5). Left: AMPT 5 < 6 < 11 fm; Middle: AMPT < 6 < 5 fm; Right: hijing 5 < fe < 11 fm. The bands correspond to 
systematic uncertainties from nonflow fit errors. 




FIG. 15: (Color online) Relative strength of average nonflow to flow {S)/{v'^{2}) as a function of the minimum »?-gap size 
(with a maximum r?-gap of 5). Left panel: AMPT 5 < 6 < 11 fm. Middle panel: AMPT < 6 < 5 fm. Right panel: hijing 
5 < & < 11 fm. Note that for hijing the denomenator ^^{2} is not flow but away-side nonflow average strength. 



may be expected from momentum conservation of dijets 
resulting in similar near- and away-side nonflow correla- 
tions at low pt- With increasing minimum ry-gap cut, the 
ratio drops. This is also expected because the near-side 
nonflow is concentrated at narrow A77 while the away-side 
nonflow is more evenly spreaded over a wide At; range. 

V. CONCLUSIONS 

In this article, we have proposed a data-driven method 
to separate particle correlations into two parts: the A?]- 
independent part and the A?7-dependent part. The for- 
mer is likely dominated by global harmonic flows due to 
the common collision geometry, and the latter can be 
identified as nonflow due to particle correlations such as 
resonance decays and near-side intra-jet correlations. 

The method uses two- and four-particle cumulants be- 
tween particles from two pseudo-rapidity {rf) bins. It ex- 
ploits the symmetry of average flow about mid-rapidity 
in a symmetric collision system. By taking the differ- 
ence of four-particle cumulants between two pairs of rj 
bins, (riaTrjis) and {r]a, —'7/3), the average flow and the Arj- 
independent flow fluctuation effect cancel, and the At;- 
dependent flow fluctuation can be deduced. Similarly by 



taking difference of the two-particle cumulants between 
two pairs of 77 bins, the difference in A?7-dependent non- 
flow correlation (plus A7y-dependent flow fluctuation ef- 
fect) is obtained from which Ary-dependent nonflow is de- 
duced. By subtracting A?7-dependent nonflow correlation 
and A77-dependent flow fluctuation, the A?7-independent 
correlation is obtained, which is dominated by flow and 
flow fluctuations. The A?7- independent nonflow contri- 
butions, such as those from back-to-back inter-jet corre- 
lations, should still remain. However, they may be neg- 
ligible in our studied px range. 

We used AMPT to illustrate our method focusing on 
medium-central collisions. We flnd the obtained A77- 
independent correlations satisfy factorization, commonly 
expected for flow. We have checked our method with 
central AMPT events, and find that the extracted flow 
and nonflow qualitatively agree with expectations of their 
centrality dependence. The obtained A77-independent 
correlation is compared to the calculated "real" flow 
from AMPT and appears to be reasonable. Our AMPT 
model studies support the validity of our decomposition 
method. 

We have also checked our method against the hijing 
event generator. The extracted AT^-independent correla- 
tions are small and comparable to the Ayy-dependent non- 
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flow correlation. This suggests that the Ar^-dependent 
nonflow and the Ayy-independent correlation may be re- 
flecting near- and away-side correlations of dijets, one of 
the main physics in hijing. The hijing results therefore 
render further support to our decomposition method. 

Our method can be readily applied to real data anal- 
ysis at LHC and RHIC. The large pseudo-rapidity ac- 
ceptance of LHC experiments is crucial. The limited 
pseudo-rapidity coverage of the RHIC detectors may pre- 
vent our method from fully extracting the Ajy-dependent 
nonflow correlation because the region of vanishing non- 
flow is likely outside the detector acceptance. Neverthe- 
less, our method should yield valuable information on 
the pseudo-rapidity dependence of flow and nonflow in 
the limited pseudo-rapidity range. 

It is important to point out that our method essen- 
tially separates Ayy-dependent and A?7-independent cor- 
relations, not particularly flow and nonflow. The sepa- 
rated correlations may be identified as flow and nonflow, 
if the A77-dependent flow fluctuation effect and the Arj- 



independent nonflow effect are small. These appear to 
be the case at low pT for the models we studied. With 
increasing px the A?7-independent nonflow effect should 
increase. Therefore, it would be also interesting to study 
our method with high-pT particles in relativistic heavy- 
ion collisions. We leave such studies to future work. 
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